Soft Matter

We probe the structural and rheological properties of ionically crosslinked alginate, a model biopolymer hydrogel, using microscopy, rheology and viscosity dependent molecular probes. This combination of techniques enables the quantification and correlation of microstructure, microviscosity, bulk viscoelasticity and yielding dynamics. By adjusting the stoichiometric ratio, R, between the alginate biopolymer and its cation crosslinks, we observe a transition from a homogeneous network-like structure to a coarse bundle-like structure at high (R>0.67) stoichiometric ratio. Intriguingly, these bundle-like structures have distinct and counter-intuitive rheological properties. Using molecular probes, we observe a continuous decrease in microviscosity that is correlated with a decrease in bulk elastic modulus and an increase in energy dissipation. This is accompanied by a transition in the hydrogel yielding under strain from a sharp, well-defined yield point to a continuous ductile-like yielding. We ascribe these surprising transitions to the looser intermolecular interaction between alginate biopolymers in the bundle-like state, as previously predicted by x-ray scattering experiments. These findings reveal new and counter-intuitive structure-property relations that demonstrate high crosslink concentration does not necessarily translate to optimal mechanical performance. SignificanceAlginate is a polysaccharide biopolymer naturally found in brown seaweed cell walls. Extracted alginate forms ionically crosslinked hydrogels that are strong, flexible and increasingly valuable in the biomedical, packaging and food industries. A large part of the utility of these hydrogels stems from the ease with which their mechanics can be tuned through adjusting the stoichiometry between alginate and its ionic crosslinks. However, little is known about how the material properties of alginate hydrogels, in particular their rheology and dynamics, are affected by microscale structural transitions at high stoichiometric ratios. Here, we use a multi-modal approach to describe and correlate hydrogel material properties and to demonstrate that increased polymer crosslinking can, counterintuitively, sometimes weaken hydrogel performance.

Successful intracellular delivery requires a nanocarrier to adhere to the cell membrane, be wrapped and internalized and finally released into the cytoplasm before it is degraded. While the initial steps have been intensively studied, the intracellular release of a nanocarrier has received comparatively less attention and remains poorly understood. To understand the physical mechanisms that govern intracellular release, we investigate the escape of nanocarriers (radii 4.5{sigma} and 7{sigma} ) from a lipid vesicle using coarse-grained simulations. We represent a nanocarrier using the metaparticle model and the membrane using the three-beads Cooke-Deserno model. Starting from a fully endocytosed state, we systematically vary the repulsion between the nanoparticle and the lipid head groups and compare three release scenarios: passive release, localized inside-out activity from a drug-like source confined within the particle and activity distributed over the entire nanoparticle. An elastic-energy analysis shows that in a tensionless membrane escape requires crossing a large pore-expansion barrier of {approx} 68kBT for the smaller particle and {approx} 110kBT for the larger, making passive release thermodynamically improbable. Our results show that the smaller particle nonetheless escapes above a threshold repulsion, whereas the larger one remains trapped in a metastable state in which local membrane pores form but fail to grow beyond the particle diameter to allow successful release. Adding activity lowers the nucleation barrier and promotes earlier pore opening, but its effectiveness depends on how it is applied, i.e., a localized internal source assists release of the smaller particle yet is insufficient to systematically free the larger one, while activity distributed across the whole nanoparticle drives complete release even in the larger case. These results show that productive escape is governed by whether the resulting pore can evolve into a geometrically accessible release pathway, providing physical design principles for intracellular delivery by nanocarriers

Mechanical properties of biological tissues, driven by passive and active forces, play a vital role in several processes ranging from development to cancer metastasis. However, the dynamical responses of cells in tissues, subject to mechanical deformations such as shear and the associated rheological properties, are not well characterized. Here, we use three-dimensional agent-based models for normal and cancer tissues to investigate their responses to simple shear as a function of cell stiffness and stochastic active forces. In the normal epithelium, with uniform strength of active force, the yield stress as a function of shear rate follows the Herschel-Bulkley form over a range of cell volume fraction. Strikingly, the shear rate dependence and the elasticity-dependent changes in the yield stress fall on master curves upon suitable scaling. To model cancer-like behavior, a certain fraction (Np) of cells was chosen to have enhanced activity and decreased stiffness. As Np increases, the extent of collective cell movement decreases, transitioning from affine (collective) to non-affine (individualistic) movement, a finding that is in accord with imaging experiments. Simulations of a model of a stiff solid tumor, with radius Rs embedded in normal tissue, show that as Rs increases, the yield stress increases. Interestingly, the cells migrate collectively as Rs increases. A Gaussian Mixture Model (GMM) and a mean field theory quantitatively account for the simulation as well as experimental results on cancerous, non-cancerous, and a mixture of these two types. The combined theoretical and experimental study establishes that heterogeneity in stiffness and activity determines non-affine movements in normal and cancer tissues.

Liquid-liquid phase separation (LLPS) driven by intrinsically disordered protein regions underlies the formation of membraneless organelles. The role of these structures in cell biology is an active area of research, yet their physical understanding remains less developed than that of protein assemblies with well-defined, rigid tertiary structures. The amino acid composition of disordered regions, particularly the abundance of charged residues, is known to modulate phase behavior. These factors govern accessible conformations, the assembly of compact sticker clusters, and the porosity of the resulting polymer networks. In contrast, the spatial distribution of side-chain volumes remains a comparatively understudied sequence determinant. Here, we employ a model that captures the patterning of side-chain volumes and interaction strengths without specifying chemical identity to show how sequence blockiness modulates the phase-nucleation propensity, connectivity, and compactness of polypeptide-like polymer assemblies. Using a library of polymers with identical composition but permuted residue arrangements, we show that clustering of bulky residues generates cooperative sticker domains and identify the threshold blockiness required to stabilize polymer networks. Coarse-grained simulations and mean-field theory reveal that this effect arises from an interplay among steric repulsion, solvent entropy, and higher-valency interactions, whereby blocky patterns promote the formation of multi-node percolated networks. Segregation of stickers and spacers is therefore one of the sequence parameters governing LLPS which may be useful for interpreting the evolution of disordered protein sequences, high-throughput proteome screening, and engineering of stimuli-responsive biomaterials. These rules may be exploited to predetermine condensate stability, viscoelasticity, and molecular exchange rates without altering chemical functionality.

Phase separating biopolymers are often described by the sticker-spacer model, which accounts for non-attractive linkers that connect discrete attractive moieties. "Associative polymer" theories describe these systems as two coexisting fluids of stickers and spacers, respectively. In such theories the stickers and spacers are decoupled in the sense that sticker binding does not affect spacer conformations or vice versa. Informed by coarse-grained simulations, we develop a phenomenological mean-field theory to determine the microscopic coupling between the sticker and spacer degrees of freedom. Accounting for linkers connecting adjacent spacers introduces free energy contributions due to conformational entropy of the spacer and the binding combinatorics of the stickers. These contributions have opposite dependencies on spacer length because longer spacers allow a greater volume to search for binding partners but incur a greater entropy cost upon bond formation. The length dependent terms from these contributions cancel leaving a length-independent elastic term arising from the fact that a small amount of spacer stretching lowers the free energy by increasing the number of possible binding partners. The cancellation of spacer length from the theory justifies the decoupling approximation and supports mean-field treatments of heterogeneous spacers.

Cell membranes are decorated by dense, polymeric macromolecules such as proteins with bulky intrinsically disordered regions, glycosylated proteins, and glycopolymers. Experiments have illustrated that the curvature of cellular membranes, which plays a role in diverse biological processes, can be regulated by the local density of these membrane-tethered polymers and their response to environmental changes. However, the mechanisms underlying the feedback between solution conditions, membrane curvature, and the local density and molecular weight of membrane-tethered polymers remain largely unexplored. In this work, using a combination of theory and experiment, we show that the ionic strength of the surrounding solution and the degree of ionization both play important roles in regulating the conformation of the membrane-tethered polymers, which in turn influences their membrane curvature generation capabilities. Specifically, we show that for negatively charged, disordered proteins and glycans, the thickness of the polymer brush decreases with ionic strength but increases with grafting density and degree of ionization, resulting in enhanced sensitivity to membrane curvature and an amplified ability to drive membrane curvature at lower ionic concentrations. Our findings reveal a mechanism by which membrane-tethered polymers generate membrane curvature to minimize free energy, with the extent of curvature generation modulated by the polymer grafting density and layer thickness, as well as by the ionic strength of the surrounding solution. The implications of these findings are broad, suggesting that perturbations to biopolymer charge, such as sialylation of glycans and phosphorylation of disordered proteins, could actively modulate membrane curvature induction in diverse cellular processes from endocytosis to cell motility.

Wound repair is an important mechanism to preserve tissue integrity in organisms after injury. However, why different tissues exhibit different mechanisms to repair wounds is a long-standing question that remains unanswered. In this work, we theoretically explore the role of the purse string, an actomyosin contractile cable used by tissues to close small wounds. Does the tissue 3D geometry influence the efficiency of the purse string in driving wound closure? Using a 3D biophysical model, we study in silico tissues with the same cell volumes but different aspect ratios, ranging from squamous to thick and tall tissues. The model predicts that taller cells are easily deformed by the purse string. In contrast, very squamous cells require a very strong purse string that might demand additional cellular mechanisms to close the gap. These findings establish a theoretical framework to predict the optimal biophysical mechanisms of wound healing in different tissues. Graphical abstractCells of different aspect ratios can be observed in a range of organisms with different function and mechanics. The wound healing efficiency of the purse string increases with the cell aspect ratio in our theoretical exploration. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/743165v1_ufig1.gif" ALT="Figure 1"> View larger version (23K): org.highwire.dtl.DTLVardef@d44ab0org.highwire.dtl.DTLVardef@1737cbaorg.highwire.dtl.DTLVardef@101b5d4org.highwire.dtl.DTLVardef@1487f26_HPS_FORMAT_FIGEXP M_FIG C_FIG

Nanoplastic interactions with cell membranes can influence particle uptake, translocation, and biological effects, making their molecular description important for assessing the environmental and biological consequences of plastic pollution. However, molecular simulations generally represent nanoplastics using monodisperse polymer chains, whereas nanoplastics formed by environmental degradation are polydisperse. Here, we construct 4 nm polyethylene terephthalate (PET), polyethylene (PE), and polystyrene (PS) nanoparticles with lognormal molecular-weight distributions reported for degraded polymers, and with melt-like chain packing and entanglement. Using Martini coarse-grained molecular dynamics, we simulate their interactions for 10 microseconds with three lipid bilayers of increasing complexity. PE and PS insert rapidly and remain largely intact, whereas PET remains surface-associated and progressively releases individual chains. Fluid, polyunsaturated lipid environments undergo greater deformation and promote substantially more chain release, up to 56% of the chains for PET and 43% for PS. Release begins with the shortest chains and progressively extends to longer chains, while semicrystalline PE remains largely intact. In contrast, a size-matched monodisperse PS particle releases no chains and minimally perturbs the membrane. These results identify the molecular-weight distribution as a key determinant of nanoplastic-membrane interactions, and indicate that uniform model particles underestimate the release of polymer chains into biological membranes.

Lipid nanoparticles (LNPs) are one of the leading platforms for delivering nucleic acid therapeutics, yet their efficacy is limited by physicochemical interactions with the extracellular matrix (ECM) that trap the particles before they reach target cells. PEGylated nanoparticles mitigate these interactions by forming a protective steric layer on their surfaces. However, there is a lack of a predictive tool that gives mechanistic insights about how PEG surface density governs the underlying interaction and results in enhanced diffusive transport of LNPs through the ECM. Here, we present a multiscale hierarchical computational framework that couples all-atom constant pH molecular dynamics (CpHMD) with a highly coarse-grained model of the complete LNP within a crosslinked hyaluronic acid (HA) network. These atomistic simulations resolve the free energy of interaction between the LNP surface and HA chains across varying PEG lipid compositions, and integrate these free energy profiles to inform the coarse-grained simulations of LNP transport through the matrix. This work highlights that even a slightly PEGylated surface depletes the near-contact shell between the LNP and HA chains, which disrupts their adhesive interactions. These protective PEG layers produce a sharp, non-linear enhancement in LNP diffusivities, with just 1% PEG increasing the diffusivity nearly eight-fold relative to bare LNPs, which remain trapped in the matrix structures. This work provides a quantitative estimate of how PEG surface density governs LNP transport through the ECM, which offers predictive guidance for engineering LNP surface properties in target-specific drug delivery.

Recent experiments show that protein condensates sitting on opposite surfaces of a flat lipid membrane move together and prefer to overlap, even though they cannot touch each other. This points to an indirect, membrane-mediated interaction. Two mechanisms could be responsible: a curvature-induced interaction, which is energetic in origin, and a fluctuation-induced interaction, which is entropic. Here we study both with coarse-grained molecular dynamics simulations, using Cookes implicit-solvent lipid model together with a generic bead-spring polymer model for the condensate. We compute the potential of mean force between two condensates across the membrane. For condensates of the same size, full overlap is unfavorable, and the pair instead settles into a partially overlapping state that bends the membrane into an S-like shape. When the two condensates differ strongly in size, full overlap becomes favorable. We explain this with a simple geometric picture. The condensate wets the membrane as a thin film and imposes curvature only along its rim, while membrane tension flattens the membrane under its interior. The resulting ring of curvature can trap a smaller condensate on the opposite side. We also compare the bending undulations and the effective bending modulus of a bare membrane, a membrane with one condensate, and a membrane with condensates on both sides. A wetting condensate suppresses the undulation modes and stiffens the membrane, but whether this makes overlap entropically favorable remains inconclusive. Our results indicate that the coupling is driven mainly by curvature, and that it depends on the wetting mechanism and on the membrane tension.

The actin cytoskeleton drives shape changes and transport in cells and supports mechanical signal transmission. However, how cells control and use active cytoskeletal flows at the mesoscale is not well understood. We reconstituted an active cytoskeleton in water-in-oil emulsion droplets of Xenopus laevis egg extract and observed, above a critical droplet size, the emergence of a 3D radially convergent steady-state flow of polymeric actin, maintained by continuous actin turnover. The flow condensed lipid-rich cellular debris into a centered inclusion. Steady-state F-actin density and flow velocity profiles roughly collapse onto scale-invariant master curves. This behavior can be explained by a physical model representing the network as an isotropic active viscous fluid with a percolation threshold. The contracting network behaves as an active swimmer with complex internal dynamics that centers itself and the central inclusion inside the droplets without physical boundary attachment. Active contraction, crosslinking and polymerization dynamics in an actin network can thus generate cell-scale flow patterns that sense the confining geometry and external signals and exert forces that are likely sufficient to move and localize organelles in cells.

Inefficient endosomal escape is a crucial barrier to intracellular delivery of nucleic acid therapeutics using lipid nanoparticles (LNPs). The use of ionizable lipids (ILs) has significantly improved cargo delivery efficiency, yet the physical basis of their interaction with endosomal membranes and their role in endosomal escape remain unclear. It has been suggested that, as ILs become cationic during endosomal acidification, electrostatic affinity promotes fusion of the LNPs with the endosome. In this paper, we propose an additional mechanism in which ILs are redistributed from LNPs to host membranes, modulating the elastic properties and curvature of the membrane, lowering the energetic threshold for endosome disruption. To test this, we quantified the spontaneous curvature of clinically relevant ILs and ATP-binding lipids and measured the membrane mechanics of giant unilamellar vesicles (GUVs) with an endosome-relevant composition at endosome-relevant pH. Small-angle X-ray scattering (SAXS) measurements reveal that the incorporation of ILs and ATP lipids into endosome-mimetic membranes shifts the spontaneous curvature towards more negative values. Micropipette aspiration experiments indicated a decrease in the apparent area compressibility modulus of membranes doped with ILs and ATP lipids. In addition, membranes showed enhanced fluctuation amplitudes and altered relaxation behavior, consistent with membrane perturbations associated with lipid insertion and pH- or ATP-driven destabilization. Under conditions promoting the partitioning of ILs or ATP-binding lipids, we further observed reduced bending rigidity and increased heterogeneity in membrane tension. Together, these results support a model in which ILs (as well as newly developed ATP-binding lipids) partition into endosomal membranes, softening the membrane and generating local curvature frustration that facilitates endosomal disruption during the natural acidification process. By quantitatively linking lipid composition with changes in membrane elasticity and fluctuation dynamics, this work provides a biophysical framework for understanding how lipid redistribution may contribute to endosomal escape and improve delivery efficiency.

Lipid membranes undergo liquid-liquid phase separation to form micron-scale domains. Several physical parameters affect this phase transition. For example, a decrease in temperature causes membranes to demix, as does an increase in hydrostatic pressure. However, measurements to determine how tension affects membrane phase separation have yielded conflicting results. Experiments that have applied osmotic pressure differences to a population of vesicles have reported an increase in the membrane's miscibility transition temperature. Conversely, experiments that have applied micropipette aspiration or substrate stretching to single membranes have reported a decrease. Here, we find that both osmotic pressure and micropipette aspiration can increase transition temperatures. We discuss how membrane pores and hidden areas in membranes present challenges to researchers seeking to quantitatively convert experimentally measured osmolarity differences into membrane tensions. We show that challenges of comparing data from different osmotic pressure experiments can be mitigated by renormalizing osmolarity differences, specifically by dividing by the exterior osmolarity. We discuss our results in the context of existing theoretical predictions and in light of four known effects of increasing tension on vesicle membranes: 1) a reduction in hidden area of tubes and aggregates, 2) a reduction in out-of-plane thermal fluctuations, 3) an increase in the area per lipid, and 4) the formation of pores. First, hidden area can explain why vesicles can sustain high osmolarity differences. Next, suppression of thermal fluctuations and increases in the area per lipid may account for shifts in miscibility transition temperatures. Finally, pores can explain time dependences. Overall, our results highlight the need for new theory and simulation that unify predictions of how transition temperatures vary over all four regimes of membrane tension.

We investigate the mechanical principles underlying fruiting body morphogenesis in Dictyostelium discoideum. Quantitative shape analysis based on the Young--Laplace law, together with AFM indentation measurements, indicate surface tension as the dominant tissue-scale force acting on the culminating fruiting body. Based on this observation, we construct a hydrodynamic phase-field model with tunable surface and interfacial tensions, and analyze its behavior numerically. Our results show that, once a stalk begins to form, the elevation of the cell mass arises naturally through a dewetting process. Through quantitative comparisons with experimental measurements, we identify the mechanical conditions required for detachment from the substrate and for establishment of the characteristic morphology of the culminating fruiting body. Together, our model analysis highlights the importance of stalk-tip elongation and tissue-scale surface and interfacial tensions in the construction of large-scale three-dimensional tissues.

Cell monolayers are active, orientationally ordered materials whose mechanics can depart from near-equilibrium behavior. Existing continuum theories often inherit assumptions from equilibrium liquid-crystal physics, motivating explicit nonequilibrium formulations. Here, we develop a minimal continuum model of two-dimensional monolayers based on odd mechanics, accounting for broken symmetries and nonequilibrium dynamics. We show that nematic order can support an odd viscous modulus generated by broken time-reversal symmetry and spatial anisotropy, without chirality. The model reproduces half-integer defect motion and stress profiles in Madin-Darby canine kidney (MDCK) monolayers, as well as defect-associated cell accumulation and depletion in neural progenitor and ovarian mesothelium systems. Finally, we estimate the viscous-moduli tensor from measured stress, velocity, and orientation fields in MDCK monolayers and identify a nonzero odd modulus. Our results show that odd mechanics provides a minimal framework for nonequilibrium cell monolayers, complementing conventional active-nematic theories.

Clathrin-mediated endocytosis (CME) is the primary pathway for internalization of extracellular and membrane cargo in eukaryotic cells. It is characterized by a patch of the plasma membrane invaginating and pinching off to generate a cytoplasmic vesicle. Recently, experimental evidence has begun to accumulate in support of proteins assembling into liquid-like condensates at nascent CME sites, but potential functional consequences are not understood. Given that condensate capillary forces can influence membrane bending, we constructed complementary analytical and numerical model systems to probe the mechanical influence of a protein droplet on endocytosis. Model evaluations and analysis reveal that a droplet can regulate the progression of endocytosis by stalling invagination until the droplet grows past a threshold size. After the droplet-mediated brake is released, the droplet provides a force that can drive the membrane bud toward scission. Thus, we propose a ``mechanical timer'' mechanism as a function for protein phase separation in cellular membrane trafficking.

Fingertip friction plays a central role in tactile exploration and object manipulation. During sliding, tangential force depends jointly on the real contact area and the interfacial shear stress, both of which can be influenced by sliding conditions. However, changes in fingertip friction are often interpreted primarily through changes in real contact area, whereas the accompanying changes in interfacial shear stress remain less well characterized. This gap is especially relevant for electrostatic surface haptic displays, which modulate fingertip friction by applying a voltage between the finger and the touch surface. Here, we experimentally quantify the mean interfacial shear stress of a sliding fingertip on an electrostatically actuated touchscreen using simultaneous measurements of tangential force and optically resolved real contact area. Ten participants performed sliding trials across three speeds and three normal forces with and without electrostatic actuation. Interfacial shear stress increased with speed and decreased with normal force; in both cases, these trends arose because real contact area varied more strongly than tangential force. Electrostatic actuation further reduced interfacial shear stress, as increasing voltage produced a larger increase in real contact area than in tangential force. These findings show that interfacial shear stress varies systematically with sliding conditions and electrostatic actuation, clarifying how changes in real contact area and interfacial shear stress combine to shape fingertip-surface friction.

Chronic wounds (principally pressure sores, venous leg ulcers and diabetic foot ulcers) are a drain on global health services and remain a major area of unmet clinical need. Chronic wounds are characterised by a bacterial biofilm (densely aggregated colonies of bacteria encased by a matrix of extracellular polymeric substances), which hinders innate immune response and can prevent wound healing. Bioactive glass (BG) fibres doped with antimicrobial metal ions, such as silver, can offer a promising treatment for chronic wound infections, where silver is well known for its antimicrobial activity against a range of pathogens and is commonly used in wound dressings. We first present a system of non-linear partial differential equations to model the treatment of a chronic wound biofilm infection with BG fibres. The BG fibres are assumed to have two mechanisms of action against the biofilm: physical disruption of the top layers of the biofilm by the BG fibres; and release of antimicrobial silver ions from the BG fibres, which then diffuse into the biofilm and can kill the bacteria. Treatment-associated parameters are estimated from in vitro experimental data using a combination of least-squares minimisation and Approximate Bayesian Computation (ABC). Sensitivity investigations are performed on other parameters that cannot currently be calculated experimentally to investigate their influence on treatment efficacy. We thus predict key parameter regimes that should lead to biofilm eradication, crucially informing the future design of metal-doped BG fibres to maximise treatment efficacy. Author summaryChronic wounds are a huge drain on global health services and will become even more problematic due to an ageing population. Current treatment methods are often unsuccessful, where treatment failure is exacerbated by the presence of a biofilm infection. Biofilms consist of communities of bacteria that adhere to the wound surface and produce extracellular polymeric substances, which can protect the bacteria by acting as both a physical and chemical barrier. More recently, there has been a focus on biofilm-based wound care, where the aim is to firstly eradicate the biofilm infection, which then enables wound healing to occur naturally. Our aim is to produce a novel treatment that can target and eradicate the biofilm infection, followed by directly assisting the wound healing. Bioactive glass (BG) fibres doped with silver offer a promising treatment as they have both anti-biofilm effects and can also stimulate the wound healing process. Here, we restrict attention to their anti-biofilm properties. By developing a mathematical model, we can predict treatment outcomes under several different scenarios, the results of which can then be utilised during design of the BG fibres. Using this combination of computational and experimental approaches, we reduce both the cost and time of optimising this promising treatment.

The rising incidence of invasive fungal diseases emphasizes the need for novel therapeutic strategies, including membrane-targeting antifungal agents, which require representative lipid models for detailed molecular-level studies. In this work, we propose a consensus quinary fungal plasma membrane model based on lipidomic literature data, specifically PC:PE:PI:PA:PS phospholipid model with ratio of 44:29:13:8:6. Using a bottom-up approach, we characterized the biophysical properties of this system - with particular emphasis on mechanical parameters such as bending rigidity and area compressibility - by combining molecular dynamics simulations with experimental flicker-noise and ATR-FTIR spectroscopies. Furthermore, we investigated the effect of two key non-phospholipid components: ergosterol and triacylglycerols. Biophysical analysis revealed that DPPI and its specific interactions with DSPS induced the most substantial deviations in baseline membrane parameters, particularly area per lipid, membrane thickness, and area compressibility, while DSPS influenced bending rigidity change and DLiPA primarily affected lipid packing defects. In addition, ergosterol and TGs were found to influence all of the investigated parameters to different degree. Notably, the overall biophysical profile of the proposed FPMM closely mimicked that of natural vesicles derived from yeast lipid extracts, establishing this model may provide a reliable platform for studying fungal membrane biophysics and lipid-targeting interactions.

Biofilm mechanics depend strongly on growth history, hydration, nutrient availability, and interfacial attachment, yet conventional rheology requires harvesting and transfer that can perturb the as-grown structure. Here, we introduce bioRheoFlow, a 3D-printed rheometer-mounted platform for cultivation and mechanical characterization without sample transfer. The platform combines evaporation compensation with static or axial nutrient flow, and benchmark measurements validate its rheometric performance. Applied to Pseudomonas putida and Pantoea sp., bioRheoFlow reveals that biofilms cultivated in situ under axial flow exhibit substantially higher elastic moduli than conventionally grown and scraped samples. Despite similar linear-viscoelastic stress ranges, the species show distinct nonlinear responses: P. putida exhibits a lower crossover stress and more abrupt failure, whereas Pantoea sp. sustains larger stresses and deformations. Creep suggests contributions from bulk deformation and biofilm-substrate interactions. bioRheoFlow provides a controlled framework for probing growth-conditioned mechanics in living biofilms.